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ABSTRACT 

Aims. To compute the chemical evolution of spiral bulges hosting Seyfert nuclei, based on updated chemical and spectro-photometrical 
evolution models for the bulge of our Galaxy, to make predictions about other quantities measured in Seyferts, and to model the 
photometric features of local bulges. The chemical evolution model contains updated and detailed calculations of the Galactic potential 
and of the feedback from the central supermassive black hole, and the spectro-photometric model covers a wide range of stellar ages 
and metallicities. 

Methods. We computed the evolution of bulges in the mass range 2 x 10' - 10" Mq by scaling the efficiency of star formation and 
the bulge scalelength as in the inverse-wind scenario for elliptical galaxies, and considering an Eddington limited accretion onto the 
central supermassive black hole. 

Results. We successfully reproduced the observed relation between the mass of the black hole and that of the host bulge. The 
observed nuclear bolometric luminosity emitted by the supermassive black hole is reproduced only at high redshift or for the most 
massive bulges; in the other cases, at z a rejuvenation mechanism is necessary. The energy provided by the black hole is in most 
cases not significant in the triggering of the galactic wind. The observed high star formation rates and metal overabundances are easily 
achieved, as well as the constancy of chemical abundances with the redshift and the bulge present-day colours. Those results are not 
affected if we vary the index of the stellar IMF from x = 0.95 to x = 1.35; a steeper IMF is instead required in order to reproduce the 
colour-magnitude relation and the present A'-band luminosity of the bulge. 

Conclusions. We show that the chemical evolution of the host bulge, with a short formation timescale of ~ 0. 1 Gyr, a rather high 
efficiency of star formation ranging from 11 to 50 Gyr"' according to the bulge mass and an IMF flatter with respect to the solar 
neighbourhood, combined with the accretion onto the black hole is sufficient to explain the main observed features of Seyfert galaxies. 

Key words, galaxies: Seyfert - galaxies: abundances - galaxies: active - galaxies: evolution 



1. Introduction The most widely accepted explanation for the luminosity 

• • emitted by an AGN, is radiatively efficient gas accretion onto 

. ^ ' The outstanding question of the co-evolution of Active Galactic a central supermassive BH. The outflows from AGNs can pro- 

k> " Nuclei (AGNs) and their host galaxies has received consider- foundly affect the evolution of the host galaxy, e.g. by quench- 

5_j able attention in the past decades, since various pieces of evi- ing or inducing the star formation (e.g., see Ciotti & Ostriker 

d dence pointed to a link between the formation of supermassive 2007, and references therein). The mutual feedback between 

black holes (BHs) and the formation and evolution of their host galaxies and QSOs was used as a key to solve the shortcom- 

spheroids: for example, the usual presence of massive dark ob- ings of the semianalytic models in galaxy evolution, e.g. the fail- 

jects at the centre of nearby spheroids (Ford et al. 1997; Ho ure to account for the surface density of high-redshift massive 

1999; Wandel 1999); the coiTelation between the BH mass and galaxies (Blain et al., 2002; Cimatti et al., 2002) and for the a- 

the stellar velocity dispersion of the host (for quiescent galax- enhancement as a function of mass (Thomas et al., 2002), since 

ies, Ferrarese & Merritt 2000; Gebhardt et al. 2000a; Tremaine it could provide a way to invert the hierarchical scenario for the 

et al. 2002; for active galaxies, Gebhardt et al. 2000b; Ferrarese assembly of galaxies and star formation (see e.g. Monaco et al., 

et al. 2001; Shields et al., 2003; Onken et al. 2004; Nelson et al. 2000; Granato et al., 2004; Scannapieco et al., 2005). 

2004) or its mass (Kormendy & Richstone 1995; Magorrian et xhe study of the chemical abundances of the QSOs was 

al. 1998; Marconi & Hunt 2003; Dunlop et al. 2003); the similar- first undertaken by Hamann & Ferland (1993), who combined 

ity between light evolution of quasar (QSO) population and the chemical evolution and spectral synthesis models to interpret the 

star formation history of galaxies (Cavaliere & Vittorini 1998; n y/C iv and N v/He ii broad emission line ratios, and found 

Haiman et al. 2004); the establishment of a good match among out that the high metallicities and the abundance ratios of the 

the optical QSO luminosity function, the luminosity function of broad-line region are consistent with the outcomes of the models 

star-forming galaxies and the mass function of dark matter halos for gi^nt elliptical galaxies (Arimoto & Yoshii, 1987; Matteucci 

(DMHs) atz ~ 3 (Haenheltet al., 1998). & Tornambe, 1987; Angeletti & Giannone, 1990), where the 

timescales of star formation and enrichment are very short and 

the initial mass function (IMF) is top-heavy. In the same year. 

Send offprint requests to: balleroOoats . inaf . it Padovani & Matteucci (1993) and Matteucci & Padovani (1993) 
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employed the chemical evolution model of Matteucci (1992) 
to model the evolution of radio-loud QSOs, which are hosted 
by massive ellipticals, following in detail the evolution of sev- 
eral chemical species in the gas. They supposed that the mass 
loss from dying stars after the galactic wind provides the fuel 
for the central BH and modeled the bolometric luminosity as 
Lboi - rjMc^, with a typical value for the efficiency of 77 = 0.1 
and were successful in obtaining the estimated QSO luminosities 
and the observed ratio of AGN to host galaxy luminosity. Then, 
they studied the evolution of the chemical composition of the gas 
lost by stars in elliptical galaxies and spiral bulges for various el- 
ements (C, N, O, Ne, Mg, Si and Fe), and found out that due to 
the high star-formation rate (SFR) of spheroids at early times the 
standard QSO emission lines were naturally explained. The rel- 
atively weak observed time dependence of the QSO abundances 
for t > I Gyr was also predicted. The model of Matteucci & 
Padovani (1993) still followed the classic wind scenario, where 
the efficiency of star formation decreases with increasing galac- 
tic mass and which was found to be inconsistent with the cor- 
relation between spheroid mass and a-enhancement (Matteucci 
1994). Moreover, Padovani & Matteucci (1993) pointed out that 
if all mass lost by stars in the host galaxy after the wind were 
accreted by the central BH, the final BH mass would be up to 
two orders of magnitude larger than observed. 

Other works (Friaga & Terlevich 1998; Romano et al. 2002; 
Granato et al., 2004), which had a more refined treatment of gas 
dynamics, limited their analysis of chemical abundances to the 
metallicity Z and the [Mg/Fe] ratio and their correlation with the 
galactic mass. 

All these studies were mainly devoted to studying the co- 
evolution of radio-loud QSOs and their host spheroids, which 
are elliptical galaxies. Now we want to extend the approach of 
Padovani & Matteucci (1993) to AGNs hosted by spiral bulges, 
with a more recent chemical evolution model for the bulge with 
the introduction of the treatment of feedback from the cen- 
tral BH and a more sophisticated dealing of the accretion rate. 
Since Seyfert nuclei are preferentially hosted by disk-dominated 
galaxies (Adams, 1977; Yee, 1983; MacKenty, 1990; Ho et al. 
1997) our study can be applied to this class of objects. 

The paper is organized as follows: in §2 we illustrate the 
chemical and photometrical evolution model, in §3 we show our 
calculations of the potential energy and of the feedback from su- 
pernovae (SNe) and from the AGN, in §4 we discuss our results 
concerning the black hole masses and luminosities, the chem- 
ical abundances and the photometry, and in §5 we draw some 
conclusions. 

2. The evolution model 

2.1. Chemical evolution 

The model on which we base our analysis is essentially the re- 
cent bulge chemical evolution model from Ballero et al. (2007a), 
which was successful in reproducing the most recent measure- 
ments of metallicity distribution (Zoccali et al., 2003; Fulbright 
et al., 2006) and evolution of abundance ratios (Origlia et al., 
2002; OrigHa & Rich, 2004; Origlia et al., 2005; Rich & Origlia, 
2005; Fulbright et al., 2007; Zoccali et al., 2007; Lecureur et al., 
2007) of the bulge giants. Here we resume its main features. 

- The gas is supposed to be well-mixed and homogeneous at 
any time (instantaneous mixing approximation). 

- The star formation rate (SFR) is parametrized as 

i//ir,t)ocvo-'' (r,t). (1) 



Mt (Mo) 


y(Gyi-') 


Re (kpc) 


taw (Gyr) 


2x IC 


11 


1 


0.31 


2 X 10'" 


20 


2 


0.27 


10" 


50 


4 


0.22 



Table 1. Features of the examined models: bulge mass (first col- 
umn), star formation efficiency (second column), bulge effective 
radius (third column). The table also reports the time of occur- 
rence of the galactic wind (fourth column). 



In this expression v is the efficiency of star formation (a pa- 
rameter whose meaning is that of the inverse of the timescale 
of star formation), o-gg, is the gas surface mass density and k 
is an index whose value can vary between 1 and 1 .5 without 
giving rise to remarkable differences. We chose ^ = 1 to re- 
cover the star formation law of spheroids. In our formalism 
the SFR is in units of MoGyr 'pc^^. 

- The bulge forms by accreting gas from the Galactic halo at 
an exponential rate: 

Mi„f cx e-"\ (2) 

where t is a suitable collapse timescale. The metallicity Zacc 
of the accreted gas is of the order of < lO^^Z©, but it is 
easily shown that the results are not significantly affected if 
we consider Z^cc - 0. 

- The Type la SN rate is computed like in Matteucci & Recchi 
(2001) following the single degenerate scenario of Nomoto 
etal. (1984). 

- Stellar lifetimes (Kodama, 1997) are taken into account in 
detail. We adopted the stellar yields of Frangois et al. (2004), 
which were constrained in order to reproduce the observed 
chemical abundances of the solar neighbourhood in the two- 
infall model of Chiappini et al. (2003). 

- The adopted IMF is constant in space and time and has the 
shape of a multi-part power law: 

(f>(M) oc M"<'+-'^'' (3) 

where the subscript refers to different mass ranges. Its nor- 
malization is performed by assuming a stellar mass range of 
0.1 -8OM0. 

In the reference model of Ballero et al. (2007a), the param- 
eters which allow to best fit the metallicity distributions and the 
[of/Fe] vs. [Fe/H] ratios measured in the bulge giants are the fol- 
lowing: V = 20 Gyr t = 0.1 Gyr, and two slopes for the IMF, 
namely xi = 0.33 for 0.1 < M/Mq < 1 (in agreement with the 
photometric measurements of Zoccali et al., 2000) and X2 - 0.95 
for 1 < M/Mq < 80. We also considered the case X2 - 1.35 
(Salpeter, 1955) for comparison. 

This model holds for a galaxy like ours with a bulge of 
Mh - 2 X IO'^Mq. We are going to predict the properties of 
Seyfert nuclei hosted by bulges of different masses, therefore 
some model parameters will have to be re-scaled. We choose 
to keep the IMF constant and to scale the effective radius and 
the star formation efficiency following the inverse-wind sce- 
nario (Matteucci, 1994). The possibility of changing the infall 
timescale with mass is not explored in the present paper In Table 
[T]are reported the adopted parameters for each bulge mass. 

2.2. The spectro-photometric model 

By matching chemical evolution models with a spectro- 
photometric code, it has been possible to reproduce the present- 
day photometric features of galaxies of various morphological 



S.K. Ballero et al.: Evolution of Seyfert galaxies 



3 



types (Calura & Matteucci 2006, Calura et al. 2007a) and to per- 
form detailed studies of the evolution of the luminous matter 
in the Universe (Calura & Matteucci 2003, Calura et al. 2004). 
In this paper, by means of the chemical evolution model plus a 
spectro-photometric code, we attempt to model the photometric 
features of galactic bulges. All the spectro-photometric calcula- 
tions are performed by means of the code developed by Jimenez 
et al. (2004); this model is based on new stellar isochrones com- 
puted by Jimenez et al. (1998) and on the stellar atmospheric 
models by Kurucz (1992). The main advantage of this photo- 
metric code is that it allows one to follow in detail the metal- 
licity evolution of the gas; this is possible thanks to the large 
number of simple stellar populations calculated by Jimenez et al. 
(2004) by means of new stellar tracks, with ages between 1 x 10*" 
and 1.4 x 10'° yr and metalhcities ranging from Z = 0.0002 to 
Z = 0.1. 

Starting from the stellar spectra, first we build simple stellar 
population (SSP) models consistent with the chemical evolution 
at any given time and weighted according to the assumed IMF. 
Then, a composite stellar population (CSP) consists of the sum 
of different SSPs formed at different times, with a luminosity at 
an age to and at a particular wavelength A given by: 



Jo Jz, 



iffito - t)LsspA^, to - t)dZdt 



(4) 



where the luminosity of the SSP can be written as: 



P,^(Z,?o-0= 4>{M)h{Z,Mjo-t)dM 



and where /i(Z, M, to - t) is the luminosity of a star of mass M, 
metallicity Z and age fo - f; Z, and Zj are the initial and final 
metalhcities, M„„„ and M,nax are the smallest and largest stellar 
mass in the population, (p{M) is the IMF and ipit) is the SFR at 
the time t. Peletier et al. (1999) have shown that dust in local 
bulges is very patchy and concentrated in the innermost regions, 
i.e. within distances of ~ \QQpc. They have also shown that dust 
extinction efl'ects are negligible at distances of ~ IRefj- In ad- 
dition, we compare our photometric predictions to observational 
results largely unafi'ected by dust, such as the ones by BalceUs & 
Peletier 1994. For these reasons, in all our spectro-photometric 
calculations, we do not take into account dust extinction. 

3. Energetics and AGN feedback 

The bulge lies in the potential well of the Galaxy, which con- 
sists of both luminous and dark matter. In the model of Ballero 
et al. (2007a) we calculated the binding energy AEg of the bulge 
gas following Bertin et al. (1992), i.e. by treating the bulge as 
a scaled-down two-component elliptical (thus ignoring the disk 
contribution). We also assumed that the thermal energy of the 
bulge interstellar medium was mainly contributed by the explo- 
sion of Type I and Type 11 supemovae. In particular, if we call 
E,h the gas thermal energy, we supposed that at the time tew, 
when the condition 



Ethitow) = ^Egitow) 



(6) 



is satisfied, a wind develops and the star formation is suppressed. 
In the reference model of Ballero et al. (2007a), this event does 
not have a great impact on the predicted chemical evolution, 
since it occurs when most of the gas has already been processed 
into stars. However, the high-metalUcity tail of the metallic- 
ity distribution was somewhat overestimated without the galac- 
tic wind. 



In the present exploration not only do we adopt a more real- 
istic disk galaxy model to better estimate the binding energy of 
the gas in the bulge, but we also consider the additional contri- 
bution to the gas thermal energy given by the BH feedback, so 
that: 



Erh(t) - E,h^SN{f) + E, 



th,AGN' 



it). 



(7) 



The various simphfying assumptions adopted in the evaluation 
of the gas binding energy and of the BH feedback will be dis- 
cussed in some detail in the following sections. 

3.1. The gas binding energy 
If we define 



and 



0.1 



0(M)Y(r, ; - TM)RM{t - TM)dM 



(8) 



(9) 



where Rm is the return mass fraction (i.e. the fraction of mass 
in a stellar generation that is ejected into the interstellar medium 
by stars of mass M; see Tinsley, 1980) and tm is the lifetime of 
a star of mass M, then, before the galactic wind, the mass of gas 
in the bulge evolves at a rate 



(5) Mg{t < taw) = Minf + M, - Y(r, t) - Mbh, 



(10) 



i.e. gas is accreted from the halo, is restituted by stellar mass 
loss and is subtracted by star formation and black hole accre- 
tion. After the wind, star formation vanishes and we suppose 
that, due to the development of a global outflow, the infall is ar- 
rested. Therefore, this equation reduces to: 



M„{t > tew) = M. - Mbh - Mn 



(11) 



where Mw is the rate of mass loss due to the galactic wind. We 
assume that all the gas present at a given time is lost, so that 
Mw = M, -Mbh- In the present treatment the galaxy mass model 
is required in order to estimate the binding energy of the gas 
in the bulge, a key ingredient in the establishment of the wind 
phase. In particular, the current galaxy model is made by three 
different components, namely: 

- A spherical Hernquist (1990) distribution representing the 
steUar component of the bulge: 



Ph(r) = 



Mh 



rb 



2n r(r + rt,)^ ' 
GMb 
r + ri, 



(12) 



where Mi, is the bulge mass, (^t is the bulge potential and rb 
is the scale radius of the bulge, related to its efi'ective radius 
hyRe ^ I. Sri,. 

A spherical isothermal dark matter halo with circular veloc- 
ity v^: 



pDMir) 



AnGrl 



ro 



(13) 



where tq is an arbitrary scale-length. 
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A razor-thin exponential disk with surface density 



R/Rj 



(14) 



where Mj is the total disk mass, R is the cylindrical radius, 
and Rd is the disk scale radius. As well known the gravita- 
tional potential of a disk in general can be expressed by us- 
ing the Hankel-Fourier transforms (e.g., Binney & Tremaine 
1987): however, as shown in Appendix|A] under the assump- 
tion of a spherical gas distribution, the contribution to the gas 
binding energy can be easily calculated without using the ex- 
plicit disk potential. 

In fact, we assume that the gas distribution before the es- 
tablishment of the galactic wind is spherically symmetric, and 
parallel to the stellar one, i.e. 

2n r(r + rb) 

In order to estimate the energy required to induce a bulge 
wind, we define a displacement radius r,, and we calculate the 
energy required to displace at r, the gas contained at r < r, (while 
maintaining spherical symmetry). We adopt a fiducial value of 
r, - 3Re; the calculated values of the binding energy do not 
change significantly for r, ranging from 2Re to IQR,,. A spheri- 
cally symmetric displacement (while not fully justified theoreti- 
cally) allows a simple evaluation of the gas binding energy, and 
it is acceptable in the present approach. Consistently with the 
assumption above, 

AEg = AEgh + AEgDM + AEgd, (16) 

where the various terms at the r.h.s. describe the gas-to-bulge, 
gas-to-DM and gas-to-disk contributions. Elementary integra- 
tions show that 



AE, 



and 



Jo 



GMgMb 



rb 



3(1 H- 5)3 



(17) 



AE. 



gDM - 



47T f pg(r)[(:>DM(n) - <iDM(r)]p-dr 
Jo 

ln(l+6)- 



where 6 = rjrb, while 



1+6 



AE 



gd - 



GMgMd 

Rrl 



xAE 



gd. 



(18) 



(19) 



and the function AE^^ is given in Appendix. Of course, AEg is 
linearly proportional to the gas mass in the bulge. 

In Fig.[T]we see that, for a IO^'Mq bulge like ours, the dom- 
inant contribution arises from the dark matter halo, whereas the 
bulge and disk contributions are comparable, both about one or- 
der of magnitude smaller than the dark matter halo one. This 
differs from previous calculations (Ballero et al. 2007a) in the 
way that the bulge contribution is reduced by almost one order 
of magnitude. The same is true for bulges of other masses. We 
must consider moreover that in the inside-out scenario for the 
Galaxy formation (Chiappini et al. 1997) the disk will probably 
form much later than the bulge, so its contribution to the poten- 
tial well during the bulge formation could be negligible. Thus, 
what we explore is the extreme hypothesis that the disk has been 
in place since the beginning. 




time (yr) 

Fig. 1. Time evolution of the different contributions to the gas 
binding energy in the bulge of a Milky Way like galaxy: dark 
matter halo (dashed line), bulge (dotted line), and exponential 
disk (soHd Hne). In particular", Mb - 2 x IO'^Mq, R^ - 2 kpc, 
Vc = 200 km/s, Md = 1O"M0 and Ra = 4.3 kpc. 



3.2. Feedback from supernovae 

The cumulative thermal energy injected by SNe is calculated as 
in Pipino et al. (2002). Namely, if we call RsNia/ii{t) the rate of 
Typela/II SN explosions: 



Eih,SN(t) - E,h^sNIa{t) + E,h^sNIl{t), 

where 



^-'th,SNIa/II 



(t) 



= r e(t-t')RsNia/n{t')dt'eTg. 
Jo 



(20) 



(21) 



The evolution with time of the thermal content e of a SN 
remnant, needed in equation above, is given by (Cox, 1972): 

s _ / 7.2 X lO^Ofo erg for 0<tsN< tc, 

'^''"^ - \ 2.2 X W^e^itsNltc)-' '^ erg for ts^ > f„ ^^^^ 

where eo is the initial blast wave energy of a SN in units of 10^' 
erg, assumed equal for all SN types, ts n is the time elapsed since 
explosion and is the metallicity-dependent cooling time of a 
SN remnant (Cioffi et al., 1988): 



1.49x10^60^/' V^V'^" 



yr. 



(23) 



In this expression ^ = Z/Zq, and no is the hydrogen number 
density. 

3.3. Accretion onto and feedback from the central black hole 

In our phenomenological treatment of BH feedback, we only 
considered radiative feedback, thus neglecting other feedback 
mechanisms such as radiation pressure and relativistic particles, 
as well as mechanical phenomena associated with jets. From 
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this point of view we are following the approach described in 
Sazonov et al. (2005), even though several aspects of the physics 
considered there (in the context of elhptical galaxy formation) 
are not taken into account. In fact, we note that these phenomena 
can be treated in the proper way only by using hydrodynamical 
simulations. 

We suppose that the bulge gas is fed into the spherically ac- 
creting BH at the Bondi rate Mb- However, the amount of ac- 
creting material cannot exceed the Eddington limit, i.e.: 



Mbh = tmn(MEdd,MB) 



(24) 



The Eddington accretion rate, i.e. the accretion rate beyond 
which radiation pressure overwhelms gravity, is given by: 



(25) 



where t] is the efficiency mass-to-energy conversion. In general, 
0.001 < T] < 0.1; we adopt the maximum value t] = 0.1 (Yu & 
Tremaine, 2002). The Eddington luminosity is given by 



Lsdd = 1.3 X 10' 



,46 Mbh 
lO^Mn 



erg s 



(26) 



The Bondi accretion rate describes the stationary flow of gas 
from large distances onto the black hole, for a given gas temper- 
ature and density (see Bondi, 1952), and is given by 



Mb = ^nRgPBCs, 
where 

GMBHjJ-mp 



Rr = 



with 



2ykT 



,^ 1 Mbh 



, (io«k) 



and pb (the gas density at Rb) can be estimated as 



Pe I Re 



(27) 



(28) 



(29) 



(30) 



In the code we adopt y = 1 (isothermal flow). If we assume that 
all the gas mass is contained within IR^, the mean gas density 
within Re is given by 



3M, 

%7lRl' 



(31) 



redshift 



6 5 4 3 



1 b- ■ 
0.1 



© 0.01 
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0.0001 



= 2 X 109 J, 
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/ 



< 
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time (yr) 
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a' 0.1 ^ 
0.01 
0.001 



Mfc = 10" Mg 



/ ' 

/ 



108 



IQi' 
time (yr) 



101 



Fig. 2. Time evolution of the Eddington (dashed Une) and Bondi 
(dot-dashed line) accretion rates for bulges of various masses. 

The thicker lines indicate the resulting accretion rate, which is 
assumed to be the minimum between the two. 



GMdMi,Wu 

■ M^vt. - 



6rb 



Rd 



(34) 



where the term Wm is given in Appendix. 

The bolometric luminosity emitted by the accreting BH is 
then calculated as: 



The equihbrium gas temperature can be estimated as the 
bulge virial temperature. 



^ f = 3.0 X 10^ K 



200 km s" 



(32) 



Lboi = nc^MBH 



(35) 



Finally, we assume that the energy released by the black hole 
is the integral of a fraction / of this luminosity over the time- 
step: 



where cr is the one-dimensional stellar velocity dispersion in the 
bulge, which is given by: 



'3Mb 



(33) 



The virial potential trace Wh for the bulge is obtained by sum- 
ming the contribution of the three galaxy components, i.e. 

Wb = Wbb + WbDM + Wbd 



-'th,BH 



pt+At 



Lboidt ^ fLboi^t 



(36) 



The value of / can vary between and 1; we assume / = 0.05 
(Di Matteo et al. 2005). 

The seed black hole mass was modifled in a range 5x10^-5 x 
10"* Mq without any appreciable change in the results. Therefore, 
we adopted a universal seed black hole mass of 10^ Mq. 
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M 
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>> 

a 
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1056 
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0.5 1 

time (Gyr) 



Fig. 3. Energy balance as a function of time for bulges of vari- 
ous masses. The figure shows the gas binding energy (solid line) 
compared to the thermal energy released by supernovae (dashed 
line) and the accreting black hole (dot-dashed line). 



4. Results 

4. 1 . Mass loss and energetics 

Fig.|2]shows the evolution of the Eddington and Bondi accretion 
rates with time and redshift. The redshift was calculated assum- 
ing a ACDM cosmology with Hq - 65 km s"' Mpc"', Qm = 0.3 
and Q.A = 0.7, and a redshift of formation of z/ ^ 10. We can see 
that the history of accretion onto the central BH can be divided 
into two phases: the first, Eddington-limited and the second, 
Bondi-limited. Most of the accretion and fueling occurs around 
the period of transition between the two phases, which approx- 
imately coincides with the occurrence of the wind although it 
extends for some time further in the most massive models. The 
details of the transition depend on the numerical treatment of the 
wind, however since the gas consumption in the bulge is very 
fast, and the Bondi rate depends on the gas density, we can ex- 
pect that the results would not change significanlty even if we 
modeled the galactic wind as a continuous wind. The BH mass 
is essentially accreted, within a factor of two, in a period ranging 
from 0.3 to 0.8 Gyrs, i.e. 2 to 6% of the bulge lifetime, which we 
assume to be 13.7 Gyr. 
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Fig. 4. Star formation rate as a function of time and redshift for 
bulges of various masses. The peak value of the Mt, = 2 x 10''Mq 
case is ~ 4OM0/yr. 



Fig.[3]compares the different contributions to the thermal en- 
ergy, i.e. the feedback from SNe and from the AGN, with the po- 
tential energy. We see that only in the case of Mt, - 2x 10"^ Mq 
does the BH feedback provide a thermal energy comparable 
to that produced by the SN explosions before the onset of the 
wincQ. Therefore, in the context of chemical and photometri- 
cal evolution, the contribution of the BH feedback is negligible 
in most cases, unless we assume an unrealistically large frac- 
tion of the BH luminosity to be transferred to the ISM. Note 
that this conclusion is also supported by hydrodynamical simu- 
lations specifically designed to study the effects of radiative BH 
feedback in elliptical galaxies (Ciotti & Ostriker 1997, 2001, 
2007; Ostriker & Ciotti 2005). Also Di Matteo et al. (2003) 
showed that it is unlikely that black hole accretion plays a cru- 
cial role in the general process of galaxy formation, unless there 
is strong energetic feedback by active QSOs (e.g. in the form of 
radio jets). 

4.2. Star formation rate 

Fig. ID shows the evolution of the global SFR in the bulge as a 
function of time and redshift. The break corresponds to the oc- 
currence of the galactic wind. It is evident that in the case of the 
most massive bulges it is easy to reach the very high SFRs of 
a few times lOOOMo/yr inferred from observations at high red- 
shifts (e.g. Maiolino et al., 2005). It is worth noting that this re- 
sults matches very well the statement of Nagao et al. (2006) that 
the absence of a significant metallicity variation up to z ^ 4.5 im- 
plies that the active star-formation epoch of QSO host galaxies 
occurred at z > 7. There is evidence for some tiny downsizing 
(star formation ceases at slightly earlier times for larger galax- 



' The Bondi accretion rate does not vanish even if the thermal en- 
ergy of the ISM is larger than the potential energy, due to the fresh gas 
provided by the stellar mass losses of the aging stars in the bulge. 
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Fig. 5. Final black hole masses as a function of the bulge 
mass predicted by our models (triangles). The short dashed- 
long dashed line represents the measured relation of McLure & 
Dunlop (2002), the dotted line the one of Marconi & Hunt (2003) 
and the dashed hne the one of Haring & Rix (2004). 



ies; see TablefTli. We stress that we are making predictions about 
single galaxies and not about the AGN population, and we only 
want to show that it is possible to achieve such high rates of star 
formation in a few Myrs. 

4.3. Black hole masses and luminosities 

In Fig.|5]we show the final BH masses resulting from the accre- 
tion as a function of the bulge mass. The predicted BH masses, 
which are reached in a few hundred Myrs, are in good agree- 
ment with measurements of BH masses inside Seyfert galaxies 
(Wandel et al., 1999; Kaspi et al., 2000; Peterson, 2003). This is 
a valuable result, given the simplistic assumptions of our model, 
since in this case we did not have to stop accretion in an artifi- 
cious way as in Padovani & Matteucci (1993); in fact, the BH 
growth at late times is limited by the available amount of gas, 
as described by Bondi accretion. The figure suggests an approx- 
imately linear relation between the bulge and BH mass, as first 
measured by Kormendy & Richstone (1995) and Magorrian et 
al. (1998). Therefore, Seyfert galaxies appear to obey the same 
relationship of quiescent galaxies and QSOs, as already stated 
observationally by Nelson et al. (2004). 

There were claims for a non-linear relation between spheroid 
and black hole mass (Laor, 2001; Wu & Han, 2001); however, 
measurements of Marconi & Hunt (2003) re-established the di- 
rect proportionality between the spheroid mass Mjp/, and Mbh- 
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also in good agreement with recent estimates from McLure & 
Dunlop (2002) and Dunlop et al. (2003): 
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Fig. 6. Evolution with time and redshift of the bolometric lumi- 
nosity emitted by the accreting BH for bulges of various masses. 
The dotted lines represent the range spanned by observations 
(see text for details). 



and from Haring & Rix (2004): 
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The three relations of McLure & Dunlop (2002), Marconi & 
Hunt (2003) and Haring & Rix (2004) are thus reported in the 
figure for comparison. The agreement between measurements 
and predictions is rather good (within a factor of two). 

In Fig.|6]we plot the predicted bolometric nuclear luminosi- 
ties Lboi for the various masses, compared with luminosity es- 
timates for the Seyfert population (see e.g. Gonzalez Delgado 
et al., 1998; Markowitz et al., 2003; Brandt & Hasinger, 2005; 
Wang & Wu, 2005; Gu et al., 2006). Of course, Uoi is pro- 
portional to the calculated mass accretion rate (Fig. |2]i. In the 
first part, the plots overlap because the Eddington luminosity de- 
pends only on the BH mass, and it is independent on the galaxy 
mass; on the contrary, the Bondi accretion rate is sensitive on 
the galaxy features. The break in the plot corresponds to the 
time when the Bondi accretion rate becomes smaller than the 
Eddington accretion rate (see Eqs.l25]andl27]l. and the break oc- 
curs later and later for more massive galaxies, which therefore 
keep accreting at the Eddington rate for longer timescales. This 
helps explaining why the outcoming final black hole mass is pro- 
portional to the adopted bulge mass. 

We see that the luminosities near the maximum are repro- 
duced by models of all masses. The same is not true for local 
(z - 0) Seyferts; only the most massive model yields a bolomet- 
ric luminosity which lies in the observed range while for other 
masses the disagreement reaches a factor of ~ 100. To simply 
shift the epoch of star formation of the less massive Seyferts, 
would require unrealistically young bulges (< 1 Gyr), whereas 
one of our main assumptions is that spiral bulges are old systems. 
Padovani & Matteucci (1993) assumed that the BH accreted the 
whole mass lost from evolved stars. If we calculate the fraction 
MbIM* at any given time, we see that its value is very close to 
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0.01. This explains why they obtained the correct nuclear bolo- 
metric luminosity of local radio-loud QSOs, but severely overes- 
timated the mass of the resulting BH, which led them to suggest 
that the accretion phase should last for a period not longer than a 
few 10'^ years, at variance with the present work. Moreover, it is 
not physically justified to assume that all the mass expelled from 
dying stars falls onto the black hole. 

One way to overcome this problem is to suppose that the 
smallest bulges undergo a rejuvenation phase, which is possible 
because spiral bulges are not isolated systems, but an interaction 
with their surrounding disks can be triggered in several ways 
(bar instabilities, minor mergers, fly-by's, and so on). Thomas 
& Davies (2006), using models combining recent star formation 
with a base old population, found that the smallest bulges must 
have experienced star formation events involving 10 - 30% of 
their total mass in the past 1-2 Gyr Such episodes may help 
explain as well the presence of the Galactic bulge Mira popula- 
tion, whose calculated age is 1 - 3 Gyr (van Loon et al., 2003; 
Groenewegen & Blommaert, 2005), and of young star and star 
clusters in the very center of the Galaxy (Figer et al., 1999; 
Genzel et al., 2003). A small fraction of intermediate-age stel- 
lar population was also detected in Seyfert 2 nuclei (Sarzi et 
al., 2007). 

Finally, we remark that by means of the present (non- 
hydrodynamical) modeling of the BH accretion rate and of the 
gas budget, a truly realistic accretion cannot be reproduced, 
whereas in a more realistic treatment the AGN feedback would 
cause the luminosity to switch on and off several times at the 
peaks of luminosity. We also note that at late times, when the ac- 
cretion is significantly sub-Eddington, a considerable reduction 
in the emitted AGN luminosity might result as a consequence of 
a possible radiatively inefficient accretion mode (e.g., Narayan 
& Yi 1994). This could reproduce the quiescence of black holes 
at the centre of present-day inactive spiral galaxies. 

4.4. Chemical abundances and metallicities 
Observations for QSOs and Seyferts 

Estimates of chemical abundance ratios in AGNs are not an easy 
task. The use of emission lines is subject to large uncertainties 
due to the dependence of the lines on several parameters which 
are difficult to quantify (e.g. column density, microtubulence, 
collisional excitation, and so on). This is a problem since most 
measurements of the [Fe/Mg] abundance ratio, which provides 
a clock for constraining the ages of QSOs and their timescales 
of enrichment, rely on the flux ratio of the Fe ii (UV bump) and 
Mg II (/12800) emission lines. The true physical origin of the 
Fe II UV complex was put in question by Verner et al. (2003), 
Baldwin et al. (2004) and Korista et al. (2004), who showed that 
the flux ratio does not scale directly with the abundance ratio 
due to the thermostatic effect of coolants. Absorption lines could 
be a better probe for QSO environments, especially narrow ab- 
sorption lines which avoid saturation and blending of important 
abundance diagnostics. However, these data require large signal- 
to-noise ratios and are therefore harder to obtain. Moreover, they 
are not free from uncertainties, concerning the shape of ionizing 
spectrum, the lack of ionization constraints, the unknown cover- 
age fraction and the exact location of the absorbers (see Hamann 
& Ferland, 1999 for an extensive review about QSO abundance 
diagnostics). 

Given these caveats, it is evident in any case that the emission 
spectrum of AGNs is particularly alike for a very large range of 
redshifts and luminosities (Osmer & Shields, 1999). The most 



probable explanation is a similarity of chemical abundances. In 
particular, the analysis of the Fe ii(UV bump)/Mg ii(/l2800) flux 
ratio in various redshift ranges (Thompson et al., 1999; Iwamuro 
et al., 2002; Freudling et al., 2003; Dietrich et al., 2003a; Baith 
et al., 2003; Maiolino et al., 2003; Iwamoto et al., 2004) was 
consistent with an [Fe/Mg] abundance ratio slightly supersolar 
and almost constant for redshifts out to z ~ 6.4; this result was 
also confirmed for z ~ 6 QSOs by very recent works (Kurk et 
al., 2007; Jiang et al., 2007). A weak trend with luminosity was 
detected (Dietrich et al., 2003a). 

Due to the time delay necessary for Fe enrichment from 
Type II SNe in a star formation history typical of elliptical galax- 
ies (see Matteucci & Recchi, 2001), this means that the sur- 
rounding stellar population must be already in place by the time 
the AGN shines and that, in a well-mixed ISM, star formation 
must have begun > 10^ years before the observed activity, i.e. at 
redshifts Zf > 10 (Hamann et al., 2002; Hamann et al., 2004). 

Many efforts were devoted to the measurement of nitro- 
gen lines, since due to its secondary naturfl the N abundance 
relative to its seed nuclei (e.g. oxygen) is a good proxy for 
metallicity. As abundance indicators, broad emission line ra- 
tios of N v/C IV and N v/He ii were analysed by Hamann 
& Ferland (1993). Hamann & Ferland (1999) also indicated 
N v/C IV and N v/O vi in narrow absorption line systems as pos- 
sible abundance indicators. Hamann et al. (2002) instead favored 
N iii]/0 III] and N v/(C iv -n O vi) as the most robust abundance 
indicators. The general conclusions drawn from these diagnos- 
tics are that the AGNs appear to be metal rich at all observed 
redshifts, with metallicities ranging from solar up to ~ 10 times 
solar Dietrich et al. (2003b) studied a sample of 70 high-redshift 
QSOs (3.5 < z < 5.0) and based on emission-line flux ratios in- 
volving C, N, O and He they estimated an average overall metal- 
licity of ~ 4-5 Zq for the emitting gas. A similar estimate (~5 
Zq) was drawn more recently by Nagao et al. (2006) who exam- 
ined 5344 spectra of high-redshift (z > 2) QSOs taken from the 
SDSS DR2; they also confirmed the detected trend in the N/He 
and N/C emission line ratios which suggests that there might be 
a correlation with luminosity, i.e. more luminous QSOs, resid- 
ing in more massive galaxies, are more metal rich. Bentz et al. 
(2004) suggested that some very nitrogen-enriched QSOs could 
be viewed at the peak of metal enrichment, e.g. near the end of 
their accretion phases, although their conclusions are not firm. 
Work carried out on intrinsic narrow absorption line systems 
(Petitjean & Srianand, 1999; Hamann et al., 2003; D'Odorico 
et al. 2004) confirmed that the observed N, C and Si abundance 
ratios are consistent with at least solar metallicities. In particular, 
D'Odorico et al. (2004) interpreted their observations as sugges- 
tive of a scenario of rapid enrichment due to a short (~ 1 Gyr) 
star formation burst. The amount of emission from dust and CO 
in high-redshift QSOs (Cox et al., 2002) corroborate the idea 
of massive amounts of star formation preceding the shining of 
the AGN. However, there is no need for exotic scenarios to ex- 
plain the production of heavy elements near QSOs (e.g. central 
star clusters, star formation inside accretion disks), since nor- 
mal chemical evolution of ellipticals is sufficient to this purpose 
(Hamann & Ferland, 1993; Matteucci & Padovani, 1993). 

The majority of conclusions concerning QSOs has been 
found to hold also for Seyfert galaxies, i.e. observations seem 
to confirm that most Seyfert galaxies are metal rich. 



- We recall that an element is secondary when it is produced starting 
from a seed non-primordial nucleus; therefore, its abundance is very 
sensitive on the metallicity of the gas where it is formed. 
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An overabundance of nitrogen of a factor ranging from about 
2 to 5 was first detected in the narrow line region of Seyferts by 
Storchi-Bergmann & Pastoriza (1989, 1990), Storchi-Bergmann 
et al. (1990) and Storchi-Bergmann (1991), and was later con- 
firmed by Schmitt et al. (1994). Further work by Storchi- 
Bergmann et al. (1996) allowed to derive the chemical composi- 
tion of the circumnuclear gas in 1 1 AGNs, and high metallicities 
were found (O ranging from solar to 2-3 times solar and N up to 
4-5 times solar). This trend was also measured in more recent 
works (Wills et al., 2000; Mathur, 2000). Fraquelli & Storchi- 
Bergmann (2003) examined the extended emission line region of 
18 Seyferts and claimed that the range in the observed [N ii]/[0 
ii] line ratios can only be reproduced by a range of oxygen abun- 
dances going from 0.5 to 3 times solar Rodriguez-Ardila et al. 
(2005), by a means of a multi-cloud model, deduced that a nitro- 
gen abundance higher than solar by a factor of at least two would 
be in agreement with the [N ii]+/[0 iii]-H line ratio observed in 
the narrow line Seyfert 1 galaxy Mrk 766. Finally, Fields et al. 
(2005a), using a simple photoionization model of the absorbing 
gas, found that the strongest absorption system of the narrow line 
Seyfert 1 (NLSl) galaxy Mrk 1044 has N/C > 4(N/C)0. 

In the circumnuclear gas of the same galaxy, using column 
density measurements of O vi, C iv, N v and H i. Fields et al. 
(2005b) claimed that the metallicity is about 5 times solar This 
is consistent with expectations from previous studies: Komossa 
& Mathur (2001), studying the influence of metallicity on the 
multi-phase equilibrium in photoionized gas, stated that in ob- 
jects with steep X-ray spectra, such as NLSls, such an equilib- 
rium is not possible if Z is not supersolar. Studying forbidden 
emission lines, Nagao et al. (2002) derived Z > 2.5Zo in NLSl, 
whereas the gas of broad line Seyferts tends to be slightly less 
metal rich. 

An overabundance of iron was suggested to explain the 
strong optical Fe ii emission in narrow line Seyferts (Collin & 
Joly, 2000) and could provide an explanation for the absorption 
features around ~ 1 keV seen in some of these galaxies (Ulrich 
et al., 1999; Turner et al., 1999) or for the strength of the FeKa 
lines (Fabian & Iwasawa, 2000). Lee et al. (1999), by constrain- 
ing the relationship between iron abundance and reflection frac- 
tion, showed that the observed strong iron line intensity in the 
Seyfert galaxy MCG-6-30-15 is explained by an iron overabun- 
dance by a factor of ~ 2 in the accretion disk. Ivanov et al. (2003) 
found values for [Fe/H] derived from the Mg 1 1 .50jum line rang- 
ing from -0.32 to +0.49, but these values were not corrected for 
dilution eff'ects from the dusty torus continuum, and therefore 
are probably underestimated. 

Metallicity and elemental abundances 

The mean values attained by the metallicity and the examined 
abundance ratios for z < 6 are resumed in Table dFl Note that all 
the [a/Fe] ratios are undersolar or solar 

Fig. |2] shows the evolution with time and redshift of the 
metallicity Z in solar units, for a standard ACDM scenario with 
Hq = 65 km s ' Mpc ', Qm = 0.3 and Qa = 0.7. 

It can be seen that solar metallicities are reached in a very 
short time, ranging from about 3 x 10^ to 10** years with de- 
creasing bulge mass (see Table |2]i. We notice that due to the 
a-enhancement typical of spheroids, solar Z is always reached 
well before solar [Fe/H] is attained, which occurs at times ~ 
(1 - 3) X 10^ years. Then, Z remains approximately constant 



^ We recall that throughout this paper we use the standard notation, 
i.e. for example [Fe/H] = log(Fe/H) - log(Fe/H)o. 



Mb (Mq) 


Z X lU 


Z X iU 


lU 


Z/Zg 


4.23 


6.11 


7.22 


[Fe/H] 


+0.83 


+0.96 


+ 1.07 


[O/H] 


+0.18 


+0.36 


+0.48 


[Mg/H] 


+0.38 


+0.55 


+0.69 


[Si/H] 


+0.79 


+0.98 


+ 1.13 


[Ca/H] 


+0.63 


+0.88 


+ 1.06 


[O/Fe] 


-0.68 


-0.60 


-0.59 


[Mg/Fe] 


-0.45 


-0.41 


-0.38 


[Si/Fe] 


-0.04 


+0.02 


+0.06 


[Ca/Fe] 


-0.19 


-0.07 


+0.00 


[N/H] 


+0.28 


+0.63 


+0.87 


[C/H] 


+0.15 


+0.24 


+0.30 


[N/C] 


+0.13 


+0.40 


+0.57 


f(Zo) (yr) 


10« 


6x 10' 


3 X 10' 



Table 2. Values assumed by the chemical abundance ratios and 
by the metallicity in solar units after the wind by bulges of dif- 
ferent masses. The table also reports the time by which Z© is 
reached in the different models (last line). 
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time (yr) 

Fig. 7. Evolution with time (bottom axis) and redshift (top axis) 
of the metalUcity in solar units Z/Zq. The dififerent curves denote 
different bulge masses as indicated in the upper left corner. 



for the contribution of Type la SNe and low-and intermediate- 
mass stars, before declining in the very late phases. The very 
high metallicities inferred from observations (e.g. Hamann et al, 
2002; Dietrich et al., 2003b, their Figures 5 and 6) are thus very 
easily achieved. More massive bulges give rise to higher metal- 
licities, which is in agreement with the statement that more lumi- 
nous AGNs are more metal rich, if we assume that more massive 
galaxies are also more luminous. This assumption is supported 
observationally by e.g. Warner et al. (2003) who found a posi- 
tive correlation between the mass of the supermassive BH and 
the metallicity derived from emission lines involving N v in 578 
AGNs spanning a large range in redshifts. 

The time dependencies of the abundances of the elements 
under study for each mass is shown in the Figs. [8] (iron) and|9] 
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Fig. 8. Evolution with time (bottom axis) and redshift (top axis) 
of the [Fe/H] abundance ratio in bulges of various masses, as 
indicated in the lower right corner. 



(a-elements, carbon and nitrogen). A fast increase in the abun- 
dances is noticeable at early times, as well as a weak decrease 
at later times, for all elements and all masses; after the galactic 
wind, i.e. for t > 0.3 Gyr (which corresponds to a redshift z - 6 
in the adopted cosmology) the abundances decrease by a factor 
smaller than 2. Such weak decrease occurs over a period of more 
than 13 Gyr. This can explain the observed constancy of the QSO 
abundances as a function of redshift (Osmer & Shields, 1999). 

The elements can be divided according to their behaviour 
after the wind: 

- O and Mg show moderate overabundances relative to solar 
They are essentially produced by massive stars (M > IOMq) 
on short timescales (< 10^ years), i.e. almost in lockstep 
with the star formation. After the galactic wind, therefore, 
their abundance cannot increase any more because star for- 
mation vanishes. They tend instead to be diluted due to the 
effect of stellar mass loss. In the Galactic bulge, at the time 
of the onset of the wind, Mg is more abundant than O, since 
although [O/H] is higher at earlier times, it declines very 
rapidly due to the dependence on metallicity of the adopted 
O yields (see Ballero et al., 2007a and references therein). 
Therefore, we expect the mean value attained by [O/H] after 
the wind to be smaller than that of [Mg/H], which is what we 
predict. In fact, the mean [O/H] for z < 6 ranges from -1-0.18 
to -hO.48 dex, i.e. oxygen is ~ 1.5 to 3 times solar, in good 
agreement with observations (Storchi-Bergmann et al., 1996; 
Fraquelli & Storchi-Bergmann, 2003); and the mean [Mg/H] 
ranges from -1-0.38 to -1-0.69 dex (i.e. 2.5 to 5 times solar), 
depending on the mass of the host bulge. Finally, we notice 
that both [O/H] and [Mg/H] reach solar values at times closer 
to those when Z - Zq for the different models, and therefore 
are better proxies for the metallicity in the bulge with respect 
to [Fe/H]. This is because in the bulge the metallicity is dom- 
inated by O. 



- Fe, Si and Ca are characterized by a bump immediately after 
the occurrence of the galactic wind, which leads to a signifi- 
cant increase in their abundance, reaching remarkably higher 
values than in the previous cases ([Fe/H]: 4-0.83 to 4-1.07, 
i.e. -7-10 times solar; [Si/H]: H-0.79 to H-1.13, i.e. -6-12 
times solar; [Ca/H]: -1-0.63 to 4-1.06, i.e. ~4-10 times solar). 
This bump is due to the combined effect of the discontinu- 
ity caused by the onset of the galactic wind at times of 0.21 
to 0.28 Gyr depending on the model mass, and of the max- 
imum of the Type la SN rate, occurring at ~ 0.2 - 0.3 Gyr 
In fact, these elements are produced also by Type la SNe (Fe 
mainly. Si and Ca in part, as they are also produced by Type 
11 SNe) which derive from progenitors with masses ranging 
from 0.8 to EMq. These elements are restored to the ISM on 
timescales ranging from ~ 30 Myr to a Hubble time. So, their 
abundances could in principle keep increasing even after the 
star formation has ceased. However, due to the adopted top- 
heavy IMF {x2 = 0.95), a relatively small fraction of low- 
and intermediate-mass stars was produced with respect to 
what we would expect if we had adopted a steeper IMF (see 
below for a comparison with the IMF of Salpeter, 1955), and 
therefore what we observe is a decrease of the Si, Fe and Ca 
abundances after the galactic wind. The estimates for Fe in 
Seyferts are lower than those calculated here, but we must 
take into account several factors which would lead to an un- 
derestimation of the Fe produced by stars (e.g. depletion by 
dust; see Calura et al. 2007b). 

- C, as well as O and Mg, is slightly overabundant (~1.3 to 
2.7 times solar). Although it is mostly produced by low- and 
intermediate-mass stars, in this same range of masses it is 
also used up to form N. In this case the bump is less pro- 
nounced and is mainly caused by the discontinuity induced 
by the galactic wind, which enhances the relative contribu- 
tion of low- and intermediate-mass stars. 

- Finally, a behaviour similar to Fe is expected from N, since 
the bulk of this element originates from stars with M < 
8M0. There are however two significant differences: first, 
due its secondary nature, its abundance increases much more 
rapidly than other elements with time (and with metallicity); 
and second, this very rapid increase hides the bump at the 
onset of the wind. We also adopt a primary N production 
in massive stars (Matteucci, 1986; see also Chiappini et al., 
2006). The effects of this choice are however only visible be- 
fore tew, i-e- when nucleosynthesis from massive stars is ac- 
tive; the decrease of the [N/H] ratio towards earlier times (i.e. 
lower metallicities) is less rapid. Then, after the wind, the 
plots with and without primary N essentially overlap. This 
was already observed in Ballero et al. (2006) and Ballero 
et al. (2007a). The mean values attained by [N/H] after the 
galactic wind ranges from 4-0.28 to 4-0.87, i.e. nitrogen is 
~2 to 8 times overabundant with respect to solar; however, 
[N/H] reaches overabundances of ~3-10 times solar at its 
peak. These results are in agreement with the estimates for 
Seyferts (Schmitt et al., 1994; Storchi-Bergmann etal., 1996; 
Rodriguez- Ardila et al., 2005). 

The difference of behaviour among the a-elements (O, Mg 
vs. Si, Ca) is even more evident in Fig. [TO] where the correla- 
tion with mass of the [a/Fe] abundance ratios after the wind is 
shown. A net separation is present: whereas [Si/Fe] and [Ca/Fe] 
are approximately solar, O and Mg are underabundant with re- 
spect to solar, O showing a more pronounced underabundance. 
This is also in agreement with the estimates of a slightly superso- 
lar [Fe/Mg] in QSOs and of its weak correlation with luminosity 
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Fig. 9. Evolution with time and redshift of the [X/H] abundance ratios for c-elements (O, Mg, Si and Ca), C and N in bulges of 
various masses. The discontinuity at ~ 2 x 10** yr in the [Si/H] and [Ca/H] and is due to the occurrence of the galactic wind which 
shortly precedes the maximum in the Type la SN rate. The thin lines in the [N/H] vs. time plot represent the results obtained by 
adopting primary nitrogen in massive stars as in Matteucci (1986). 
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(Dietrich et al., 2003a) Since both [Fe/H] and [Mg/H] are con- 
stant within a factor of 2 up to high redshifts (z - 6) and follow 
the same declining trend with time, we expect this relation to 
hold for most of the bulge lifetime, which is what is observed for 
QSOs (Thompson et al., 1999; Iwamuro et al., 2002; Freudling 
et al., 2003; Dietrich et al., 2003a; Barth et al., 2003; MaioHno et 
al., 2003; Iwamoto et al., 2004). Here, we predict [Fe/Mg] val- 
ues ranging from -1-0.38 to -1-0.45, the highest values correspond- 
ing to less massive galaxies where star formation stops later and 
Type la SNe have more time to pollute the bulge ISM. 

Fig. [TT] shows the variation with mass of the [N/C] abun- 
dance ratio. The figure illustrates well that this ratio is remark- 
ably sensitive to the galaxy mass, which, as we have seen in Fig. 
|2l is correlated with the galaxy metallicity. This happens because 
of the secondary nature of N, which is produced at the expense 
of C in a fashion proportional to the metallicity of the ISM. The 
mean [N/C] ratio for z < 6 ranges from H-0.13 to -1-0.57 depend- 
ing on the bulge mass, i.e. the N/C ratio is ~ 1.4 to 3.7 times 
solar. If we consider the amount of variation of the two abun- 
dances, these values are consistent with the estimates of Fields 
et al. (2005a). 

Finally, we make a brief comparison of the results obtained 
with the top-heavy IMF (x2 = 0.95) with those we obtain if 
instead we adopt X2 = 1.35 (Salpeter, 1955). In Ballero et al. 
(2007a) this IMF was excluded on the basis of the stellar metal- 
licity distribution, since with our model the Salpeter index for 
massive stars gives rise to a distribution which is too metal-poor 
with respect to the observed ones (Zoccali et al., 2003; Fulbright 
et al., 2006), and Ballero et al. (2007b) further seemed to con- 
firm this point. However, Pipino et al. (2007) showed that a hy- 
drodynamical model for ellipticals can be adapted to the galactic 
bulge and, with a spherical mass distribution, it reproduces the 
above mentioned metallicity distributions with a Salpeter IMF. 
Furthermore, as we shall show (see ^4.5l l the (B - I) colors and 
bulge A'-band luminosity are better reproduced by a Salpeter 
IMF above IMq. Therefore a brief comparison is useful. 



In Fig. [12] we show the evolution with time and redshift of 
the abundance ratios relative to hydrogen for some of the ele- 
ments considered, namely Fe, O and N. We find similar trends 
for all the elements. However, in the Salpeter case, the abun- 
dances at the wind are lower than with a top-heavy IMF, because 
the enrichment from massive stars is lower Then, those elements 
which are produced mainly by low- and intermediate-mass stars 
(in this case, Fe and N) keep increasing after the wind until they 
reach a maximum value which is generally slightly higher than 
that obtained with the top-heavy IMF, then their abundance stops 
increasing. The bump in the [Fe/H] is ~ 0.1 dex larger because 
the production of Type la SN progenitors is favored in this case. 
On the contrary, oxygen abundance decreases steadily after the 
wind because it is essentially produced by massive stars and its 
value remains below the one calculated with the top-heavy IMF. 
The abundance ratios are still almost constant for z < 6; thek 
mean values are respectively: 

- [Fe/H] = 0.79 to 1.08, i.e. 6 to 12 times solar; 

- [O/H] = -0.09 to 0.27, i.e. 0.8 to 2 times solar; 

- [N/H] = 0.29 to 0.84, i.e. 2 to 7 times solar. 

The values for Fe are thus again overestimated with respect to 
observations, but as we said, the latter are not corrected for Fe 
depletion by dust or dilution by dusty torus continuum (Ivanov 
et al. 2003). The total metallicity Z ranges from about 3.5 to 6.5 
solar, i.e. slightly less than in the top-heavy IMF case, but still in 
agreement with observational estimates, as well as the N and O 
abundances. 

We point out that the adoption of a Salpeter IMF would result 
into little change in the quantities calculated previously, i.e. mass 
accretion rate, luminosity, energetics and final black hole mass. 
The accretion rates of Eddington and Bondi do not depend on 
the adopted IMF. A steeper IMF only slightly shifts the time of 
occurrence of the wind of a few tens Myrs ahead, thus prolong- 
ing the Eddington-limited accretion phase which, as we said, is 
the phase when most of accretion and shining occurs. This will 
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Fig. 12. Evolution with time and redshift of the [X/H] abundance 
ratios for Fe, Mg and N in bulges of various masses and the adop- 
tion of a Salpeter (1955) IMF above IMq. The discontinuity at 
~ 3 X 10^ yr in the [Fe/H] and [Si/H] and is due to the occurrence 
of the galactic wind which shortly precedes the maximum in the 
Type la SN rate. 
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Fig. 13. Evolution of the predicted (U - B) (upper panel) and 
(B - K) (lower panel) colours for bulge models of three different 
masses (dotted line: Mh - 2 x lO^M©; solid line: Mb = 2 x 
10'°Mo; dashed line: Mb = IO^Mq), assuming an IMF with 
X2 = 0.95. 

lead to a larger BH mass; however, the final BH masses increase 
by only about 10%. 

4.5. Bulge colours and the colour-magnitude relation 

In this section, we present our results for the spectro-photometric 
evolution of the three bulge models studied in this paper. The 
photometrical evolution of Seyfert galaxies, which requires 
modeling of the AGN continuum and, for Type 2 Seyferts, of 
the dusty torus, is not treated at present and may be the subject 
of a forthcoming paper. 

In Fig. [13] we show the predicted time evolution of the 
(U - B) and (B - K) colours for the three bulge models studied 
in this work. In Figure 1, we assume an IMF with xi - 0.33 for 
stars with masses M in the range 0.1 < M/Mq < 1 and X2 - 0.95 
for M > IMq. At all times, higher bulge masses coiTespond to 
redder colours, owing both to a higher metallicity and to an older 
age. This is consistent with the popular "downsizing" picture of 
galaxy evolution, according to which the most massive galaxies 
have evolved faster than the less massive ones (Matteucci 1994, 
Calura et al. 2007c). 

In Fig. [141 we show how the assumption of two different 
IMFs affects the predicted evolution of the (U - B) and (B - K) 
colours, for a bulge of mass Mb - 2 y. lO'^M©. We compare 
the colours calculated with the IMF presented in ^2.11 with the 
ones calculated by assuming an IMF with xi - 0.33 for stars 
with masses M in the range 0.1 < M/Mq < 1 and X2 - 1.35 for 
M > IMq. It is interesting to note how, at most of the times, the 
assumption of a flatter IMF for stars with masses m > IMq im- 
plies redder colours. This is due primarily to a metallicity effect, 
since a flatter IMF implies a larger number of massive stars and 
a larger fraction of O (the element dominating the metal content) 
restored into the ISM at all times. 
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Fig. 14. Evolution of the predicted (U - B) (upper panel) and 
(B-K) (lower panel) colours for a bulge model of mass M/, = 2x 
lO'^M© assuming two different IMFs. The solid lines represent 
the colours computed by assuming X2 = 0.95, whereas the dotted 
lines are computed assuming X2 = 1.35. 



In Figure [15] we show the colour magnitude relation (CMR) 
predicted for our bulge models and compared to the data from 
Itoh & Ichikawa (1998, panel a), and Peletier & Balcells (1996, 
panels b and c). 

Itoh & Ichikawa measured the colours of 9 bulges in a fan- 
shaped aperture opened along the minor axis, in order to mini- 
mize the effects of dust extinction. In panel a of Fig.[T5]we show 
the linear regression to the data observed by Itoh & Ichikawa 
(1998), indicated by the solid line, whereas the two dotted lines 
express the dispersion. 

Peletier & Balcells (1996) determined the colours for a sam- 
ple of local bulges, for which they estimated that the effect of 
dust reddening is negligible. Of the sample studied by Peletier 
& Balcells (1996), here we consider a subsample of 17 bulges, 
for which Balcells & Peletier (1994) have published the absolute 
R magnitudes. This allows us to plot an observational CMR. 

The predicted CMR has been calculated by adopting two dif- 
ferent IMFs, i.e. with X2 - 0.95 and X2 = 1.35, represented in 
Fig. [15] by the solid circles and the open squares, respectively. 
From the analysis of the / vs. (B - 1) plot, we note that the adop- 
tion of a flatter IMF leads to an overestimation of the predicted 
(B - I) colours; in fact, the predictions for the flatter IMF lie 
above the upper dotted lines, representing the upper limits for 
flie data observed by Itoh & Ichikawa (1998). On the other hand, 
the predictions computed with an IMF with X2 - 1 .35 are con- 
sistent with the available observations. 

From the analysis of the R vs. (U - R) diagram, we note 
that our predictions computed assuming X2 = 1.35 are consis- 
tend with the observations, in particular for the bulge models 
of masses 2 x lO'^M© and lO'^M©. For the lowest mass bulge, 
corresponding to the highest absolute R magnitude, the (U - R) 
colour seems to be overestimated. In Table |3] we present our re- 
sults for the predicted present day colours for the bulge model 
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Fig. 15. Predicted and observed colour-magnitude relation for 
bulges. The solid line and the dotted lines are the regression and 
the dispersion of CMR observed by Itoh & Ichikawa (1998), re- 
spectively. The solid circles and open squares are our predictions 
for the three bulge models studied in this paper, computed by as- 
suming for the IMF X2 = 0.95 and X2 = 1.35, respectively. Panel 
b: R vs. (U-R) diagram. The open circles are the data by Peletier 
& Balcells (1996). Solid circles and open squares as in panel a. 
Panel c: R vs. (R - K) diagram. Open circles, solid circles and 
open squares as in panel a. 



of mass 2 x 10 "'M©, computed for two different IMFs, com- 
pared to observational values for local bulges derived by var- 
ious authors. The present-day values are computed at 13 Gyr. 
Note that, for each colour, the observational values represent the 
lowest and highest observed values reported by the authors. The 
values we predict are compatible with the observations of local 
bulges, which show a considerable spread. However, the sam- 
ple of bulges considered here seems to be skewed towards high 
mass bulges, hence it is difficult to infer the trend of the R vs. 
iU - R) for magnitudes R < -16 mag. On the other hand, the 
predictions computed by assuming X2 - 0.95 seem to produce 
too much high (U - R) colours with respect to the observations. 

Finally, the R vs. (R - K) diagram does not allow us to draw 
any firm conclusion on the slope of the IMF. 

From the combined study of the / vs. (B-I), R vs. (U-R) and 
R vs. {R - K) diagrams, we conclude that the observational data 
for local bulges seem to disfavour IMF flatter than X2 - 1 .35 in 
the stellar mass range 1 < M/Mq < 80. In Table [3] we present 
our results for the predicted present day colours for the three 
bulge models presented in this paper, computed for two different 
IMFs. The present-day values are computed at 13.7 Gyr. In the 
same table, we show the observational data of Peletier & Balcells 
(1996) and with the one of Galaz et al. (2006), who estimated 
that for their sample, dust effects should be small, with an upper 
limit on the extinction of 0.3 mag. Note that, for each colour, the 
observational values represent the lowest and highest observed 
values reported by the authors. 
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Photometric properties of tlie Galactic Bulge 





(B-R) 


iU-R) 


{R-K) 


Observed (a) 

ib) 


1.27 - 1.6 
0.5-2.18 


1.24-2.16 


2.39 - 2.97 


Predicted {x = 1.35) 
(x = 0.95) 


1.64-1.71 
1.75 - 1.83 


1.94 - 2.09 
2.20 - 2.40 


2.22 - 2.38 
2.49 - 2.68 



Table 3. Predicted colours for the Galactic Bulge assuming two 
different IMFs, compared to observational values of local bulges 
from various sources {a: Balcells & Peletier 1994; b: Galaz et al. 
2006.) 



For the model with mass comparable to the one of the 
Bulge of the Milky Way galaxy, i.e. the one with a total mass 
M = 2 X 1O'°M0, with the assumption of an IMF with X2 = 0.95 
we predict a present /sT-band luminosity Lbui,k = 4.5 x lO^L©. 
Existing observational estimates of the /iT-band luminosity of the 
bulge indicate values 0.96 x 10'" < LBui,K,obs < 1.8 x IO'^Lq 
(Dwek et al. 1995, Launhardt et al. 2002), i.e. al least a factor of 
2 higher than the estimate obtained with our model. By adopting 
an IMF with X2 = l .35, we obtain L'gJ^ f = 7.6 x IO^Lq, in better 
agreement with the observed range of values. 

5. Summary and conclusions 

We made use of a self-consistent model of galactic evolution 
which reproduces the main observational features of the Galactic 
bulge (Ballero et al., 2007a) to study the fueling and the lumi- 
nous output of the central supermassive BH in spiral bulges, as 
fed by the steUar mass loss and cosmological infall, at a rate 
given by the minimum between the Eddington and Bondi accre- 
tion rates. A realistic galaxy model was adopted to estimate the 
gas binding energy in the bulge, and the combined effect of AGN 
and supernova feedback was taken into account as contributing 
to the thermal energy of the interstellar medium. We also investi- 
gated the chemical composition of the gas restored by stars to the 
interstellar medium. Assuming that the gas emitting the broad 
and narrow lines observed in Seyfert spectra is well mixed with 
the bulge interstellar medium, we have made specific predictions 
regarding Seyfert metallicities and abundances of several chem- 
ical species (Fe, O, Mg, Si, Ca, C and N), and their redshift evo- 
lution. Finally, we calculated the evolution of the {B - K) and 
(U - B) colours, the present-day bulge colours and /T-band lu- 
minosity and the color-magnitude relation, and discussed their 
dependence on the adopted IMF. 

Our main results can be sunmiarized as follows: 

1 . The AGN goes through a first phase of Eddington-limited ac- 
cretion and a second phase of Bondi-hmited accretion. Most 
of the fueling (and of the shining) of the AGN occurs at the 
transition between the two phases, which coincides approx- 
imately with the occurrence of the wind. Within a factor of 
two, the final black hole mass is reached in a fraction of the 
bulge lifetime ranging from 2 to 5%. 

2. The peak bolometric luminosities predicted for AGNs resid- 
ing in the bulge of spiral galaxies are in good agreement with 
those observed in Seyfert galaxies (^ 10"*^ - lO'*'* erg s"'). At 
late times, the model nuclear luminosities (produced by ac- 
cretion of the mass return from the passively aging steUar 
populations) are < 10^^ erg s"\ and would be further re- 
duced by a factor ~ 0.01 when considering ADAF (or its 
variants) accretion regimes. To recover an agreement with 



the luminosities of local Seyferts, it is necessary to assume a 
rejuvenation event in the past 1-2 Gyrs. 

3. The feedback from the central AGN is not in any case the 
main responsible for triggering the galactic wind, since its 
contribution to the thermal energy is at most comparable to 
that of the supernova feedback. 

4. The proportionality between the mass of the host bulge and 
that of the central black hole is very well reproduced without 
the need to switch off the accretion ad hoc. We derive an 
approximate relation of Mbh ~ 0.0009Mfc, consistent with 
recent estimates (McLure & Dunlop 2002; Marconi & Hunt 
2003; Haring & Rix 2004). 

5. The massive amounts of star formation (several thousands 
solar masses a year) necessary to produce the observed line 
strengths, dust content and metalUcities of Seyfert galaxies 
are easily attained. Due to such high star formation rate, 
solar metallicity is reached in less than 10^ yr, and solar 
abundances for the elements we considered are reached in 
a few hundred million years. This naturally explains the high 
metalUcities and abundances inferred for Seyfert galaxies in 
analogy with QSOs. 

6. After the first ~300 million years, the interstellar medium in 
the bulge reaches overabundances of up to one order of mag- 
nitude for N, Fe, Si, Ca, up to 5 times solar for N and 3 times 
solar for C and O. This is roughly consistent with the esti- 
mates for Fe, N and O in the broad and intrinsic narrow line 
regions of Seyfert galaxies. The slightly supersolar [Fe/Mg] 
and its weak dependence on the galaxy luminosity (mass) are 
recovered. These results remain essentially unchanged if we 
adopt a Salpeter (1955) stellar initial mass function. 

7. The mean value assumed by the [N/C] abundance ratio is 
very sensitive on the galaxy mass, due to the sensitivity of 
N to metallicity. The calculated values are consistent with 
recent estimates (Fields et al., 2005a). 

8. Higher bulge masses imply redder colours, in agreement 
with the downsizing picture of galaxy evolution. A flatter 
IMF gives rise to a redder bulge. While the present-day bulge 
colours are well fitted both by a top-heavy (x2 = 0.95) and 
a Salpeter (xi = 1.35) IMF, the latter is required in order to 
achieve a better agreement with the colour-magnitude rela- 
tion and the present-day /T-band luminosity of the bulge. 

Acknowledgements. SKB wishes to acknowledge Antonio Pipino for interesting 
suggestions and discussions. 

Appendix A: The disk gravitational energy 
contributions 

In this paper, in order lo estimate the binding energy of the gas 
in the bulge, we assumed spherical symmetry for the gas distri- 
bution itself. Thus, the contribute of the spherical steUar bulge 
and the spherical dark matter halo are given by the elementary 
expressions (17) and (18). 

A more complicate case is represented by the disk-gas inter- 
action. In fact, one could suppose that a 2-dimensional integral, 
involving special functions must be evaluated, due to the disk ge- 
ometry. However, this is not the case: in fact, if each gas element 
is displaced radially from r to r,, then 

AEgd = [ pg{r)r^dr f t/Q[<*d(erO) - 0d(err)] = 

Jo ^ J-in 

= An {' Pg{r)r\Mrt) - hir)W (A. 1) 

Jo 
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where 
Mr) = 



r)dQ. 



(A.2) 



is the angular averaged potential, and is the radial unit vector. 
It is trivial to prove by direct evaluation that the angular mean of 
the potential of a generic density distribution is also the (spher- 
ical) potential of the angular averaged density. Of course, this 
result can also be obtained in a more elegant way by considering 
the angular mean of the expansion of the potential in spherical 
harmonics (e.g., Binney & Tremaine 1987), so that all the terms, 
except the monopole, vanish. An identical argument holds also 
in the case of the computation of the virial trace. In fact, from 
X = it follows that 



Wbd 



< X, V(pd > d X - -An 



d(l>d(r) 3 

r 

dr 



Pdr (A3) 



The surface density of a razor-thin axisymmetric disk in spheri- 
cal coordinates is 



S(r, <p) 
so that 
Pdir) 



Ed(rsinl?)5(cosi^) 



=- r 

47r j4; 



Ur,'&,^)da=^^. 

4n ^r 



For the exponential disk in equation (14) we then obtain 
Ur) = -'^{X-e-'^'), 

so that integrations in (Al) and (A3) reduce to 
■ h){\ - e-^l^) 



(A.4) 



(A.5) 



(A.6) 



Kir ^ + (1 



and 



1 _ e^'^l+h) / I 
'h 1? '\h 



■r 



(A.7) 



(A.8) 



e~'r^dt is the 



where h = Rd/rh,6 = r,/rh, and Ei(-x) 
Exponential Integral. 
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